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We linearize extended ADM-gravity around the flat torus, and use the associated Fock vacuum to con- 
struct a state that could play the role of a "free" vacuum in loop quantum gravity. The state we obtain 
is an element of the gauge-invariant kinematic Hilbert space and restricted to a cutoff graph, as a natural 
I/"") ' consequence of the momentum cutoff of the original Fock state. It has the form of a Gaussian superpo- 

sition of spin networks. We show that the peak of the Gaussian lies at weave-like states and derive a 
relation between the coloring of the weaves and the cutoff scale. Our analysis indicates that the peak 
weaves become independent of the cutoff length when the latter is much smaller than the Planck length. 
' By the same method, we also construct multiple-graviton states. We discuss the possible use of these 

, states for deriving a perturbation series in loop quantum gravity. 

on 
o 

o 

"q ! 1 Introduction 

cr 

! In the research on canonical loop quantum gravity (LQG) and spinfoam models 1 , one of the major open 
' problems is the development of a reliable semiclassical analysis. At present, we do not know if these theories 
contain semiclassical states which correctly reproduce the observed Einstein gravity. Nor is a perturbative 
expansion available that would allow one to calculate the scattering of low-energy excitations of such states. 
Several lines of research have led to proposals for vacuum states, tentative results on perturbations, and 
modified dispersion relations for matter: among them are approaches based on the Kodama state [5]-[10], 
spin network invariants [11, 12], linearized gravity [13, 14], coherent states [15]-[23], weaves [24]- [28] and 
general boundaries [29]. While the kinematics of the gravitational sector and the dynamics of the matter 
sector are relatively well understood, the gravitational dynamics and its semiclassical limit remain largely 
unclear. Consequently, we have no conclusive evidence that loop quantum gravity has a physically realistic 
semiclassical and low-energy limit 2 . 

The work of this article is motivated by the following question: can one obtain a semiclassical pertur- 
bation series for LQG and, if yes, how? At first sight, it may seem odd to ask this question about a theory 
that has often been characterized as the non-perturbative alternative to perturbative approaches in quantum 
gravity. It is the very failure of traditional perturbation theory that led to the loop approach to quantum 
gravity, and one of its strongest points is the fact that it does not rely on any approximative scheme for 
its definition. On the other hand, we know that theories can have both perturbative and non-perturbative 
regimes, depending on which scale we look at them. In QED, for example, perturbative expansions work 
fine at low-energies where the fine-structure constant is small, while at higher energies the coupling running 
grows and perturbation theory breaks down. It is renormalization that relates the different regimes of the 



1 For an introduction, see e.g. [1, 2] and [3, 4]. 
2 For an overview on the literature, see sec. II. 3, 



same theory, and allows us, in principle, to compute low-energy actions from actions at more fundamental 
scales. It gives us, in particular, the relation of measurable couplings at accessible scales to bare couplings 
at cutoff scales. 

When the transition from one scale to another involves strong coupling regimes, the renormalization 
procedure may require non-perturbative techniques. Thus, it can happen that a field theory is perturbative 
at low energies, while a non-perturbative renormalization is needed to compute its low-energy action. This 
suggests that the failure of perturbative gravity may not come from the coupling expansion itself, but from 
the perturbative renormalization that fails to provide us with unique couplings in the effective action. Such a 
viewpoint is supported by the work of Reuter and others (see e.g. [30]- [33]) who investigate non-perturbative 
renormalization group flows of gravitational actions. 

From that perspective, it is conceivable that loop quantum gravity is a non-perturbative theory at its 
fundamental scale, and nevertheless accessible to perturbative treatments at lower scales. While possible 
in principle, this scenario remains rather elusive, since so far the loop formalism does not dispose of the 
techniques to implement renormalization and perturbation theory. Tentative ideas on renormalization have 
been formulated by Markopoulou and Oeckl [34, 35], and there are first attempts to develop a perturbation 
theory around the Kodama state [9, 10]. Starodubtsev constructs string- like excitations of this state [7], 
and Smolin has analyzed string perturbations [36] of causal spin networks [37, 38]. 

Here, we will take a more conservative strategy and try to proceed in close analogy with ordinary 
quantum field theory. There, the working principle of perturbation theory can be roughly described as 
follows: we choose a classical background and consider only states which are semiclassically peaked around 
it. As a result, the Hamiltonian is dominated by the lowest order terms in the fluctuation — the free part 
— whereas higher orders can be treated as small corrections. The free Hamiltonian defines a linear system 
and provides a first approximation to the dynamics of the perturbative regime. Inspite of its simplicity, it 
is far from clear how this scheme should be transferred to the framework of loop quantum gravity: how 
can we generalize field-theoretic concepts to a theory, where the basic degrees of freedom are not fields, but 
labelled abstract graphs? What is the analogue of a field fluctuation in a space of networks? What tells us 
how to separate an operator on graphs into a "free" and "interaction" part? 

In the present paper, we approach this problem by constructing candidate states for a free vacuum and 
free gravitons in loop quantum gravity. We analyze these states in the hope to gain information on how 
semiclassical properties manifest themselves in the loop framework and how this could be exploited to do 
perturbation theory. 

Since we do not know how to linearize LQG itself, our approach is indirect: we know how to linearize 
ADM gravity, and we know that field and loop-like degrees of freedom are physically related — after all 
they arise from a quantization of the same classical theory. We employ this relation to translate the free 
vacuum and free gravitons of ADM gravity into states of the loop representation. In contrast to earlier 
attempts in this direction[13, 14], we arrive at states in the Hilbert space of the full non-linearized theory. 
An important feature of the construction lies in the fact that it starts from momentum-regularized states 
and translates this property into a cutoff graph of the final loop states. Thus, they can be viewed as more 
or less coarse-grained states, depending on the value of the cutoff parameter. 

Both the vacuum and gravitons take the form of Gaussian superpositions of spin networks whose graphs 
lie on the cutoff graph. We analyze the maximum of the Gaussian and find that spin networks at the peak 
have similar properties as weaves [24]: we determine their graphs and establish a relation between mean spin 
and cutoff scale. Our analysis indicates that in the limit of very small cutoff length, the peak spin networks 
become independent of the cutoff graph, have spin label 1/2 and graphs of a length scale close to the Planck 
length. That is, the graphs of these spin networks maintain an effective Planck scale discreteness, while the 
mesh of the cutoff graph becomes infinitely fine. 

The paper is organized as follows: in section 2, we linearize extended ADM gravity on a flat torus, 
quantize the reduced system and implement the linearized transformation to Ashtekar-Barbero variables. 
In this way, we obtain a free vacuum that is a functional of reduced triads or connections. Section 3 describes 
in several steps how we adapt this state to the degrees of freedom of loop quantum gravity. In section 4, 
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we use the same procedure to define loop analogues of many-graviton states. Section 5 describes the peak 
analysis. In the final section, the construction and properties of the states are summarized, and we discuss 
the similarities and differences to other proposals for vacuum states. We also mention ideas on a genuine 
loop quantization of the free vacuum which would allow one to remove the cutoff graph similarly as in the 
definition of Hamiltonian constraint and area operator. In the last subsection, we return to the issue of 
perturbation theory: we discuss how our states might be used for extracting a free part of the Hamiltonian 
constraint, and how they could be applied in renormalization. Interestingly, we arrive at an ansatz that is 
closely related to Smolin's string perturbations [36]. 



Notation and conventions 



spatial indices: a, b, c . . . = 1, 2, 3 
internal indices: i, j, k . . . = 1, 2, 3 
3- metric: g ah 

determinant of 3-metric: g 

Planck length: l p = \/fiK 

gravitational coupling constant: k = 8nG/c 3 



LQG: loop quantum gravity 

S: gauge-invariant spin network state 

S: gauge-variant spin network state 

H: kinematic Hilbert space of LQG 

Ho : gauge-invariant kinematic Hilbert space 

7~tdiff' gauge- and 3d-diff-invariant 
kinematic Hilbert space 



We use units in which c 



2 Vacuum of linearized Ashtekar-Barbero gravity 

Loop quantum gravity is based on a quantization of the so-called Ashtekar-Barbero variables — classical 
phase space variables that arise from a canonical transformation of the standard (extended) ADM variables. 
In order to arrive at some linearized version of LQG, we have, loosely speaking, six possibilities, correspond- 
ing to the different orders in which quantization (Q), linearization (L) and canonical transformation (C) 
can be applied. Let us abbreviate them by 

1. LCQ 4. QLC 

2. CLQ 5. QCL 

3. LQC 6.CQL, 

where the order of operations goes from left to right, starting with the classical ADM-theory and ending up 
with a linearized form of LQG. Combinations 4. and 5. are merely hypothetical, since the ADM-variables 
have never been rigorously quantized. Ideally, what we want is number 6., a linearization of full LQG. So 
far, however, we do not know how to do this (or, for that matter, what linearization should exactly mean 
in that case), since the degrees of freedom of the theory are quite different from fields. Hence, before that 
problem is resolved, we have to content ourselves with options 1. to 3. Actually, 1. is identical to 2. (what 
other meaning should be given to LC than CL?), so there remain possibilities CLQ and LQC: the former 
has been explored by Ashtekar, Rovelli and Smolin [13] for imaginary, and by Varadarajan [14] for real 
Immirzi parameter, when they applied a loop quantization to the linearized Ashtekar-Barbero variables. 
The route we follow in this paper takes the third variant LQC as its point of departure: we linearize 
the classical extended ADM- gravity, apply a Schrodinger quantization to it, and implement the linearized 
canonical transformation within the quantum theory. Thus, we arrive at a vacuum state which, at first 
glance, has little to do with loop-like degrees of freedom. The problem of relating this state to LQG will be 
the subject of section 3. 
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2.1 Linearization of classical extended ADM formulation 



In classical theory, linearization rests, similarly as in quantum theory, on the idea that we choose a back- 
ground, restrict attention to small deviations from it, and exploit this to give a lowest order approximation 
for the dynamics. In the following, we go through the classcial linearization of extended ADM gravity, but 
we will not motivate or derive each step. The procedure is similar to that for standard ADM and complex 
Ashtekar gravity, which has been described in the literature 3 [40, 13]. 

To keep things as simple as possible, we choose space to be the 3-torus T 3 and linearize around a flat 
background on T 3 x R. The linearization consists of the following steps: we linearize the classical constraints, 
use them to obtain the reduced phase space, and determine the Poisson brackets and Hamiltonian on it. 
Once we have obtained the reduced classical system, we will quantize it using a Schrbdinger representation 
(sec. 2.2) and, finally, apply the linearized form of the canonical transformation (sec. 2.3). 

In the extended ADM formulation (see [1]), the phase space variables are the 1-density triad 

d 

E i (x) = Ef(x)—, * = 1,2,3, (1) 
and the canonically conjugate one-forms 

K\x) = Ki(x)dx a , . = 1,2,3, (2) 

They are relatted to the 3-metric g a b and extrinsic curvature K a b by 

g a b = \E\E a E l b , (3) 
K k a = K ab E bk /^, (4) 

where E denotes the determinant det(E? ). The Poisson brackets read 

{E?(x),E b k (y)} = {Ki(x),K*(y)} = 0, (5) 

{Et(x),K k (y)} = ^S^S(x-y). (6) 
The constraints are first-class and consist of the gauge, vector and Hamiltonian constraint: 

Gjk = K a[j E^ (7) 
V a = -D h \KiE)-b h a K{E^ (8) 

C = ^(KlKl-KiKj)E*E!-yJ\E\R(E). (9) 
R(E) stands for the 3d Riemann tensor when written as a function of the densitized triad. 



Linearization around flat torus 

Choose a flat classical background triad on T 3 such that the torus corresponds to a cube with macroscopic 
side length L and periodic boundary conditions. Moreover, choose, once and for all, a coordinate system in 
which the background Ei- and K l -&e\ds read 

E a clk = 5 a k , K k a = K clab E b J/^ = 0. (10) 

We introduce the relative variables 

4 5 a k , K k = K k a -0 

3 For a detailed exposition of linearization in the Hamiltonian context, we refer the reader to [39]. 
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and adher from now on to the convention that 

e ka = e a k , K ka ^K k a . (11) 

That is, the spatial index of e can be freely moved between upper and lower right position, while the spatial 
index of K may be either in the lower or upper right position. The Poisson brackets of the relative variables 
are 

{e ab (x), K cd (y)} = {E b a (x) - S b , K c d (y)} = \ S c J d b S(x - y) . (12) 
By keeping only linear terms in e and K, we arrive at the linearized constraints 

Q ab = K [ab]^ ( 13 ) 

V a = -d b K ba + d a K, (14) 
C = 2d a 8 b e ab , (15) 

which are again first-class. 

For the phase space reduction, it is convenient to change to Fourier space. On the 3-torus, we use the 
following conventions for Fourier series: 

/(*) = "T^E 6 **/^ and (16) 

k 

f(k) = -Ljd'xe-^fix), (17) 

where the wavevector k takes values in 27r/LZ 3 . The delta functions on position and Fourier space have 
the respective transforms 

%-x') = ^E ei - fe "- ) - ( 18 ) 

k 

and 

^' = 7/ d ^ e " i( -"- ) -- ( 19 ) 

The Poisson bracket becomes 

{e ab (k),K cd *(k')} = ^S c a 5 d b 6 kig . (20) 

By imposing the linearized constraints and choosing suitable gauge conditions, we require that e ab and K ab 
are symmetric, transverse and have a constant trace. This defines our reduced phase space. We denote the 
reduced variables by e™ b and K? b d , and write the Poisson bracket again as { , }. 

When Fourier-transformed, the reduced variables can be decomposed into six zero-mode components, 
and two components for each nonzero k, corresponding to the two polarizations of gravitational waves: 

Oz) = ^(j:em^am + Y.i2e-- e ^a b (k)) , (21) 

V V \i=l k>0i=l J 

K b M = ^(E^(0ha 6 (0) + EE e --^(^^(Al)] • (22) 

V V \i=l fc>0i=l / 

More specifically, for each nonzero pair {k, —k} we choose a right-handed coordinate system s.t. one of the 
vectors, say k, points in the positive 3-direction. Then, we define the polarization tensors by 

eiab(k) ■= {5\ahb + foa$lb) , (23) 

£2ab{k) := -j=(5i a S lb - 5 2 a^2b) , (24) 

eiab(-k) := -eiab(k), (25) 

£2ab{-k) ■ = ^2ab{k)- (26) 
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It follows that 
and 



e* ab(k)tj ab(k) — §ij 



(27) 



e* ah (k) = e iab (-k). (28) 
For k = 0, we take £j a fo(0), i = 1, . . . , 6, to be an orthonormal basis in the space of symmetric 2-tensors, i.e. 

e* ab (0)e jab (0) = 6 ij . (29) 

The projection of an arbitrary two-tensor T ab onto its reduced part is given by 

(30) 



(PT) ab (x) := 4rE e--Pab cd (k)T cd (k), 
vv . 



where 



Pab cd (k) = e iab (k)ef*(k). 



(31) 



Recall that the Poisson brackets of the reduced phase space are the pull-back of the Poisson brackets on 
the full phase space. In our notation, 



This implies that 



Equivalently, we have 



{ea b d (x),K c d d (y)} = {(Pe) ab (x),(PKy d (y)}. 

{elt(k),K^(k')} = {(Pe) ab (k),(PKT d *(k>)} 

= Pa b ef (k)P cd gh (k'){e ef (k), K9 h *(k')} 

= Pa b ef (k)P cd 9 h^)^ h A,kJ 

= lP ab cd {k)h,e ■ 

{em,K*(k!)} = ^6 kA , 



(32) 



(33) 
(34) 



for polarization and zero mode components. 

Let us come to the linearized dynamics: the Hamiltonian is given by 



H = - f d 3 x N d (x)C quadr (x) 



(35) 



where C qua d r is the quadratic part of the Hamiltonian constraint (9) when evaluated on the reduced phase 
space, and N c i is the lapse density associated to the background: i.e. the lapse for which / d 3 x N(x)C(x) 
generates a flow that leaves the phase space point of the background fixed. In the case of the flat background, 
this is just N c i = ^fg. A straightforward calculation yields 



-!■ 

K J 



H = — I d x 

K 



Kt(x)KZ{x) + d c e- b \x)d c e^(x) 



'■ah 



When expressed in terms of Fourier or polarization components, the Hamiltonian reads 

H = -E[^ b :(k)K^M + k 2 eT(k)e2 d (k) 
k 

= -Y.Wtik^i^ + k'em^k) . 



(36) 

(37) 
(38) 
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The polarization components for k / describe the spatial change in e™£(x): they oscillate in harmonic 
potentials and always stay near ei(k) = 0. The zero modes are the constant part of the e^ d (x) field and 
move in a flat potential. This means that, to linear approximation, the overall shape of the torus behaves 
like a free particle. Unless the initial momentum is zero, the size of ej(0) will grow, so that at some point 
the linear approximation breaks down. This instability is due to the compactness of the torus. On R 3 , the 
zero modes are absent and the linearization stable. 



2.2 Reduced phase space quantization 

We quantize the reduced system by introducing operators e^ d (/c) and K^ d (k), and replace the Poisson 
brackets (34) and (33) by commutation relations 



ei{k),K]{g) 



= ih-P ab cd (k)5 kA , , 



The Hamilton operator becomes 

k 

= -Y,[Khk)Mk) + k 2 4(k)e t (k) 

K, L 



(39) 
(40) 

(41) 
(42) 



To represent these operators, we use a Schrodinger representation in terms of functionals of ei(k), where 
ii(k) and K l {k) act as multiplicative and derivative operators respectively: 



ei(k) -> ei(k) , Ki(k) 



-in 



d 



2de*(k) 



(43) 



Note that the coefficients e,(A;) have to satisfy the reality condition e*(k) = et{—k). That is, ej(0) is real, 
and for k / only half of the coefficients (say those for k 1 > 0) can be taken as independent variables. 
Thus, we choose the functional measure as 



/"• = {A/ 



dei(0) 



f 2 1 

n nn/ 



de ir (k) 



(44) 



The additional index r = 0, 1, denotes the real and imaginary part respectively. 

Next we specify a free vacuum for the system. We require of it that it is time-independent and peaked 
around ej(fe) = 0. The peakedness is necesssary to ensure consistence with linearization 4 . At first, one 
might think that we look for the ground state of the Hamiltonian (42): 



*G [<*(£)] = Afexp 



1 



^-E ke*(k)a(k) 



(45) 



For nonzero k ^ 0, this functional is a Gaussian around ei(k) = and satisfies our requirements. A 
dependence on the zero modes is missing, however, so *$>g has an infinite spread in ej(0), which is inconsistent 
with linearization. To remedy this, we add a Gaussian factor 



cxp 



1 



co e 2 t (0) 



UJQ > 0, 



(46) 



4 Of course, the peak property is just a minimum requirement: even then, the peak could be too wide, so that linearization 
is not applicable. 
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for each zero mode. This gives us the new state 



where 



*[ei(fc)] = AAexp 



u{k) 



k, k > 0. 
UJQ , k = . 



(47) 



(48) 



Note that under evolution by If, the spreading of the zero mode wavefunction proceeds on a time scale 
r ~ 1/uq. By choosing ljq of the order l/L or smaller, we can make the state \I/ practically time-independent 
for all microscopic processes. 

In the following, we use this state as the free vacuum of ADM gravity on T 3 x R. (Throughout the text 
we write normalization factors unspecifically as J\f and do not keep track of their precise value.) 



2.3 Canonical transformation 

The classical Ashtekar-Barbero variables are obtained by the transformation 

where (3 is the Immirzi parameter and Y l a [E] denotes the spin connection as a function of E: 



1 



K[E] = \e l3k E\ [d b Ei - d a E{ + E*E l a d b E; 



1 

+ t£ 



ijk 



a E b E 



Here, we take (3 to be real. The transformation leads to the new Poisson brackets 

{E?(x),Ai(y)} = {Et(x),ri[E}(y) + (3Ki(y)} = ^ 5 b a 5j 5{x - y) , 



or in terms of Fourier modes 



The linearization of (49) induces a canonical transformation on the reduced variables: 



At d = e a cdd c e^ + (3K, 

One may check that Af b d is again symmetric, transverse and of constant trace. 
In the quantum theory, we introduce the new operator 

At = e acd d c eX + f3K' d . 

Using that 

i£acdk c £ldb = ke 2a b an d i£acdkc^2db = keiab , 

we see that (54) is equivalent to 

Mk) = ke 2 (k) + PK^k) , 



ah 



A 2 (k) = fcei(fc) + f3K 2 {k) , 

Mo) = (3km. 



k^O. 



(49) 
(50) 

(51) 
(52) 

(53) 



(54) 
(55) 
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In the Schrodinger representation, Ai acts as 



Mk) = ke 2 (k)-ih K(3 8 



MO) 



fcei(fc) — ih 

k(3 d 
—lTi- 



2 de\{k) ' 

k(3 8 £^°> 
2 de* 2 {k) ' 



2 aef(O) • 

Up to an A-dependent phase (which we choose to be zero), eigenstates of A have the form 

2i 



<PA[ei(k)] = exp 



Hk(3 



(56) 



Within the quantum theory, the canonical transformation (53) is implemented by a unitary map 

O e i/[e]/fto e -i/[e]/R ) 
that turns the Ai-operator into a pure functional derivative in e*, i.e. 



2 0e?(fc) 



(57) 



We see from (56) that the required factor is 



e^ h :=e^\^^ke\{k)e 2 {k) \ . 



The transformed vacuum state reads 
y[ei(k)} = AAexp 
In terms of reduced Fourier components, it takes the form 
*[e$(k)]=Mexp 



2i 



Y.[^)e*{k)e l {k) - - ke\{k)e 2 {k) 



(58) 



(59) 



By doing a Gaussian integration, we can transform (58) to the ^-representation: 



*[Ai(k)]=/fexp{- — 



3/4 



2i 



(60) 



The Schrodinger representation we have defined so far is formal, since arbitrarily high momenta appear. We 
can resolve this either by a rigorous continuum formulation (employing Gaussian measures on tempered dis- 
tributions), or by using a regularization. Here, we choose the second possibility, as the existence of infinitely 
high momenta is physically questionable anyhow. We introduce an ultraviolet cutoff on the momenta, which 
we denote by A. We will sometimes refer to this as the regularized Fock representation, since it is the natural 
home for Fock-like excitations, as opposed to the polymer-like excitations in the loop representation. (Of 
course, strictly speaking we are dealing with a Schrodinger and not a Fock representation.) 
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With this adjustment, the canonically transformed and regularized vacuum becomes 



"if" E ( w ® e T(k)e- b d (k) + \ e-tik) e acd k c et{k) 



(61) 



In the regularized scheme, we consider the position space field e^ b d (x) as a function 



e^{x)-.= ^=Y,^-^{k). (62) 

k<A 



of the Fourier modes e r ^ b d (k). With that in mind, we can write the state functional also as 



, (63) 



where the kernel W\ is defined by 

(T'\ ( r. a \ — 



W A (x,y) = ^J2e i ^-y)u(k). (64) 

fc<A 

3 Transition to degrees of freedom of LQG 

At this point, we have a well-defined expression for a vacuum functional and we would like to translate it 
into a state of the LQG Hilbert space. In doing so, we want to preserve the physical properties of the Fock 
state as far as possible, which, includes, in particular, the cutoff on the momenta. The problem is that our 
state functional involves degrees of freedom that appear quite different from those of loop quantum gravity: 
on the one side, we have fields or their Fourier transforms, and on the other side abstract networks with 
spin labellings. In order to achieve a meaningful transition from Fock to loop state, it will be essential to 
find the right way to relate these degrees of freedom. 

The Fock state was obtained by a reduced phase space quantization: we imposed the linearized gauge- 
and diff-constraint on the classical level, and then quantized the reduced degrees of freedom. Thus, we 
arrived at a state which is a functional of the reduced connection (see eqn. (60)). 

LQG, on the other hand, is based on a quantization of the full phase space variables, yielding the 
kinematic Hilbert space H. The full non-linear gauge- and diff-constraint are imposed subsequently to 
give the gauge- and diff-invariant Hilbert space Ho and Hdig respectively. The configuration space, from 
which the Schrodinger representation on H is built, consists of generalized connections — distributional 
connections with support on graphs, of which ordinary connections are only a special case. States in Hdis 
are functionals of gauge- and diff-equivalence classes of such generalized connections 5 . 

Logically, we can divide this difference between degrees of freedom into three steps: 

1. reduced connection — ► connection, 

2. connection — > generalized connection, 

3. generalized connection — ► gauge- and diff-equivalence class of generalized connections. 

Our strategy for bridging this gap: We modify the Fock state ^ such that it becomes a functional 
of connections (sec. 3.1 and 3.2); then we switch from a pure momentum regularization to a combined 
momentum / triangulation based regularization which gives us f as a functional of generalized connections 
(sec. 3.3). In the final step, an averaging over the gauge- and diffeomorphism group has to be applied in 
order to arrive at a state in HdiS- In this paper, we only do the gauge-averaging explicitly, which provides 
a state in Ho (sec. 3.4). The diff-projection remains to be done. 

5 The standard procedure is to define states in Ho as gauge-invariant functionals of generalized connections, and states in 
Tidiff are constructed as equivalence classes of states in Ho under diffeomorphisms. We obtain a mathematically equivalent 
formulation by defining everything in terms of equivalence classes of connections. Then, an element of Hdm is a functional of 
equivalence classes of connections w.r.t. gauge- and diff- transformations. 
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3.1 From reduced to full configuration space 

We take the Fourier coefficients Ef(k) as coordinates for the full configuration space. The position space 
field 

Et(x) = ^=J2e i ^Ef(k) (65) 
fc<A 

should be understood as a function of the Ef(k). Recall that ei a without index rod denotes the difference 

e ia (x) = e°(x) = Et(x)-5? (66) 

between the background triad and Ef . We introduce a Schrodinger representation for functionals on the 
full configuration space: the measure is defined by 



3 3 1 



/oo 

fc<-A,/c^U i=io=ir=u ~°° 



(*), 



(67) 



and we represent the canonical commutation relations 



(cf. (52)) by setting 



E?{k)=E?(k) and ti a {k) = ° 



2 dEf*{k) ' 

Analogously, there is an ^4-representation, where Ef(k) and A l a (k) are represented by 

Ai(k) = Ai(k) and Ef(k) = ih^ ° 



2 dA*(k) 

To start with, we simplify the state (63) by dropping the /3-dependent phase factor: 



*Kl?ai)]=AAexp 



-|d 3 x| d 3 y W A (x,y)e- b d (2L)ea b d (y) 



(68) 
(69) 

(70) 

(71) 



Thus, we avoid overly long formulas in the computations that follow. The treatment with phase factor is 
discussed in section 3.5. 

Our aim is to extend the functional (71) to the full configuration space. The most simple possibility 
would be to use the projection map 

e2 d (k) = P^(k)e cd (k), (72) 
and define the extended state by the pull-back, i.e. 



*ext :=Afexp 



1 



d 3 x / d 3 y W A (x,y) (Pe) ab (x)(Pe) ab (y) 



(73) 



The problem with this state is that it has a very degenerate peak. In the vicinity of the background triad, 
this is ok because it corresponds to gauge-, diff- and time reparametrization invariance. If one goes farther 
away from the background triad, however, the linearized transformations are no longer symmetries of the 
theory. That means that if we follow long enough along the degenerate direction, we will arrive at triads 
that are very diff- and gauge- inequivalent to the chosen background, but they are still in the peak of the 
state functional due to the projector. That is a very unphysical property. 

As an alternative, we could drop the projectors in (73) and define the state as 



*ext [£?(£)] :=AAexp 



(74) 
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It implies that we make the Gaussian peak non-degenerate and throw out diff- and gauge-symmetry com- 
pletely. After the state has been transferred to the LQG Hilbert space, the lost symmetries need to be 
restored. We can do this by applying an averaging over the gauge and diffeomorphism group 6 . If we did 
the averaging in (74), it would effectively reintroduce the projection on the symmetric transverse part for 
small e, while for large e, it would establish the correct non-linear gauge- and diff-symmetry for the state. 
The problem: at the level of the loop state, gauge-averaging turns out to be rather complicated, if we start 
from (74). 

In the present paper, we use a variation of this approach: to simplify the gauge-averaging, we replace 
the triad fluctuations e ab {x) in (74) by the fluctuation of the densitized inverse metric 7 g ab (x): 



e ab (x) - \ (g ab (2L) ~ S ab ) = \ (E?(x)E b {x) - S ab ) . 



(75) 



The triad fields are 1-densitites, so g ab (x) has density weight 2. Since Ef contains only modes up to k = A, 
we can write g ab (x) also in a smeared form 



~gf := J dV E?(x)5 A (x - x')E b (x') , 
where the smearing is done with the regularized delta function 

The new state ^ ex t is defined as 
*cxt[£?(A0] := AAexp 



(76) 



(77) 



fc<A 



^-Jd 3 xJ d 3 yW A (x,y) (^j dV Ef (x)6 A (x - x')E b (x') - <5' 

J dV E?{y)5 A {y_-yjE b y)-5' 



(78) 



Note that this state is almost gauge- invariant, but not completely, due to the smearing at the cutoff scale. 
What is the justification for changing from the state (74) to (78)? With our substitution (75), the term 



becomes 



eab(x)e ab (y) 

J dV E-(x)5 A (x - x!)E b (x') - 5 ab ^j \ (J d 3 y' E a k (y)5 A (y - y^)E b k {yJ - 5 

\ (E?(x)E b {x) - 8 ab ) (Et{y)E b k {y) - 5 ab ) 

(5 + <(*)) (S b + e b (x)) - S ab ] {(5 a k + et(y)) (s b + e b k (y)) - 5 ab 

6? e b (x) + ef (x) 6 b + ef (x) e b (x)} [6 a k e b k (y) + e%{y) 6 b + e a k (y) e b (y) 

\(e b a (x)+et(x)) (e b a (y) + e a b (y)) + o(e 3 ) 
e(ab)(x) e (ab) (y) + o(e 3 ) . 



(79) 



(80) 



(81) 



6 By dropping the projector in (74), we also lost the linearized scalar constraint: it is related to the dynamics and should 
reappear when the dynamics is linearized at the level of loops. For ideas in this direction, see section 6.4. 
7 We thank C. Rovelli for suggesting this modification. 
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The basic premise of linearization is that for a suitable coupling parameter, the values of e can be divided 
into "large" and "small" fluctuations with the following property: the "large" e are irrelevant because for 
these values the state functionals under consideration are exponentially damped. The remaining "small" 
fluctuations lie near the peak of the states and are small enough that to first approximation, higher orders 
in e can be neglected relative to the leading order terms. Thus, for small e, and within the precision of 
the linear approximation, (80) and (79) are equal except for the symmetrizers — a degeneracy which is 
due to the gauge-invariance in (78). Since we intend to apply a gauge-averaging anyhow, we can ignore 
this difference. For large e, both state functionals are exponentially damped, so that again the difference 
between (80) and (79) is not important. 



3.2 "Complexifier" form 

By a Gaussian integration, we transform the state (78) to the ^-representation: 
*ext[4(*0] = M JDE exp ( ]T At(k)El{k) 



k<A 



x exp 



- / d 3 x J d 3 y W A (x,y) (/ dV E?(x)5 A (x - x')E^(x') - 8 

' j 'dV Et{y)5 K {y-y[)Ety)-8 



Using that 



j DE exp 



2i 



n n n n / d ^>® ^ e ^g^oo 

V, fci> i=la=lr=0 J ~°° / V k<A 



K k<A, fc!>0 

is the delta functional on the connection, and that the operator E^ acts like 

Kp d 



ih- 



2 dAl(k) ' 



we can write the entire expression (82) as 

1 



*«tK(fc)] = Afexp 



AflK 



d 3 x I d' 6 y W A {x,y) ( / d'^x' Ef(x)8 A (x-x')El(x') - 8 ab 



x ( / dV E?(y)8 A (y-y')E>>(y>)-8 ab 



(82) 



where 



Et{x) = n =Y.^--Et(k). 

k <A 



8(A), 
(83) 

(84) 



This form of the state is similar, but not identical, to Thiemann's general complexifier form for coherent 
states [20]. Thiemann writes coherent states as 



exp (-C 1 ) 6(A - A) 



(85) 



cl 
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where A^j is the so-called complexified connection and contains the background triad. What we do here is 
somewhat different because we leave the background triad outside the delta-function. For that reason, it is 
a slight abuse of terminology if we call our analogue of C, i.e. 

C := J d 3 xj d 3 y W A (x,y) (/ dV Ef(x)S A (x - x')^(x') - 5 ab ^j 

x (| d 3 y' E?(y)5 A (y - y')4V) - 5 ab ^j , (86) 

again a complexifier. 

3.3 Prom momentum cutoff to triangulation 

At this point, we have a well-defined expression for a state functional of connection Fourier modes. In the 
next step, we turn this into a functional of generalized connections. We do so by a change of regularization 
scheme: similarly as one changes from a momentum to a lattice regularization in ordinary QFT, we trade 
the UV cutoff on the connection for a triangulation of space. Let T A denote the simplicial complex of this 
triangulation. The connection is replaced by a generalized connection A on the dual complex 8 T^: a map 
that sends every edge e of into a group element g e = A(e). In making this transition, we want to alter 
the physical properties of the state as little as possible. The triangulation-based regularization should be 
such that it mimics the effects of the UV cutoff. For that reason, we choose T A to be regular in the following 
sense: 

When measured against the background metric 5 a b, the edges e of the dual complex should be straight 
and have lengths l e in the range 

(l-e)/ A < l e < (l + e)/ A , 

where 

l A := I (87) 

is the length scale corresponding to the cutoff A, and e is some small fixed number. 

The state we have so far consists of two parts: a delta functional of the Fourier coefficients of the 
connection and an operator acting on it. Let us first consider the delta functional: we replace it by the 
delta functional 5t*(A) on which is equal to a sum over all gauge- variant spin networks S restricted to 

5 T *(A) = J2 S(0)S*(A) (88) 

5CT A * 

In the operator part in (83), we have to replace the smeared operator product by an operator on functionals 
of A, or equivalently, by an operator on spin network states. In other words: after having quantized 

~gf = J dV E?(x)5 A (x - ^Efix') (89) 

in a Fock space manner, we will now quantize it along the lines of loop quantum gravity. 

Loop quantization of smeared inverse densitized metric 

There are several consistent ways in which one could quantize expression (89) on a triangulation. Of all the 
possibilities we will choose one that is very simple and gauge-symmetric. 

8 Given 7a, the dual complex is denned in the standard way, using the metric information of the background. 
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Before dealing directly with (89), let us give a definition for Ef(x). Spin network states are built from 
representation matrices U 3 {g e ) associated to edges e. To define Ef(x) on Uj(g e ), we think of the latter as 
the path-ordered exponential of a connection A l a (x), i.e. 

Cj)\ 



U 3 \g e ) = Ve^\^-Jds e a ( S )Al(e(s)) J 

~ f 7 V !./! (?) 

= E (^J Jd Sl jds 2 --- J d Sn e 01 ( ai )^i(e( ai ))---e a »(5 n )^,(e( an )) J n • • • J n , 



and view the triad operator as the functional derivative ihn0/2 5/5A l a (x) . With this prescription, we obtain 
that 

E?(x)Uj(g e ) = ~*y jf ds e a (s) 5(x - e(s)) U 3 (g ei ^) J t^c*)) , (90) 

where e\{x) and e2(x) are the edges which result from splitting the edge e at the point x. (When x is not 
on e, the definition of e\{x) and e 2 {x) is irrelevant, since the delta function gives zero.) The problem with 
(90) is that the original holonomy is split into two holonomies. Thus, when applied to a spin network on 
T£, the result will be a spin network which lives on the dual of a refined triangulation. If we want to stick 
to our original intention of defining the state onT A *, we have to modify the action (90) such that it leaves 
the space of spin networks on invariant. 
The choice we take is 

E?{x)U 3 (g e ) = (jT 1 ds e a (s)) (s(x - e(l)) j] U 3 (g e ) + S(x - e(0)) U 3 (g e ) J^ . (91) 

This means that -Ef-operators can only create J's at the beginnings and ends of edges. In the obvious way, 
equation (91) generalizes to an action of Ef(x) on an entire spin network state on 



V 

e of v 




E?(x)S = Y: E KM • (92) 



The sum ranges over all vertices of the spin network, and the diagrams indicate where the J's are inserted. 
The vertex- edge form factor F£ e stands for 

f^M ■= "x Of ds e{s) ) s & - ^ ■ (93) 

Given this definition of Ef(x), what would be a meaningful way to implement 

If = J dV Ef(x)S A (x - x')E^(x') ? (94) 

Roughly speaking, the smearing function S\(x — x') requires x and x' to be closer than I a- Since the cutoff 
length l\ is also the length scale of the dual edges, we translate this into the condition that, if Ef(x') inserts 
a J at a node v, Ef(x) can only insert J's at the same node. 

This still allows for the possibility that Ef(x) and E'-(x') insert J's on different edges of the same node. 
The action of the operator (94) becomes 

J dV Ef(x)S A (x - x')E^(x') S 

J dV F« ei (x)S A (x - x>)F* e2 (x!) ¥ £1 . (95) 



V edges 

ei,e2 of v 
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Due to the contraction of i-indices, this action is gauge- invariant. For a generic basis of spin networks, it is 
not diagonal, however, and that makes its use complicated. In the following, we take a simpler choice, and 
adopt the point of view that Ef(x)Ef(x') should not involve contractions of Jj's from different edges. If we 
drop these cross-terms, the action of (94) is diagonal and simply reads 



where gf b (x) is given by 



~gf(x) 



J dV Ef(x)5 A (x - x')E b (x') S = gf{x) S , 
= EE fd 3 x'F« e (x)S A (x-x!)FZ !e (x')h 2 j e (j e + l) 

v e of v 

= E E <e(x)je(je + l), 



(96) 



(97) 



v e of v 



and 



F v ab M := Jd 3 x'h 2 FZM5A(x-x!)F^ e (x') 



f 1 ds e a (s)^j Qf * ds' eV)) J dV 5{x - x v ) 5 A (x - x') 6{x' - x v ) 



= ( l -?pj £ds£ds' e a (s)e\s')5 A (0)8(x-x v ). 



(98) 



We interpret g a ^ as a densitized inverse metric that is associated to the spin network S. This metric is 
distributional and has only support on vertices of the spin network graph. Thus, we can also write 

~9i(x)=J29§\v)t(x-2Lv) (99) 



where 



9§( v ) ■= E JeUe + l) ['ds ['ds' e a (s)e b (s') 



(100) 



With the definition (96) and the replacement of the delta functional by (88), we obtain the state 



* = ^(0)exp 



SCTT 



= AT ^(O) exp 

ScT* 

In the second line, we abbreviated 



I d3x I d3y Wk{ ~ ~ & " 6<lb ) (?s® " 6ab ) 

-^-J d3x J d "y ^a(x - y) hf(x) hf(y) 



hf := gf - 5 ab . 



S* (101) 
(102) 

(103) 



Since spin networks on T£ are naturally mapped to spin networks in TC, the state (102) trivially extends to 
a state in 7i. We denote this extended state by the same formula. Note that the restriction to spin networks 
on 7^ in Ti is not ad hoc, but designed to preserve the momentum cutoff of the original Fock state. At the 
same time, the momentum regularization is not replaced completely, as it is still present in the kernel W\ 
and the factor <5a(0). 
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3.4 Gauge projection 

It remains to make the transition from (102) to a gauge-invariant state in Hq. Gauge-averaging simply 
yields 

1 



= M J2 S(0) exp 



ScT* 



4Hk 



d*x / d 6 y W A (x - y) hf(x) hf(y) 



where the sum ranges over all gauge- invariant spin networks S on T£. We introduce the coefficient 

1 



*o(S) := 5(0) exp 
and write this more compactly as 



AflK 



d 6 x d i yW k {x-y)hf(x)hf(y) 



^ =M *o(S)S*. 



(104) 

(105) 
(106) 



SCT* 



One can think of ^>o{S) as the wavefuntion of in the 5-representation. 

When expressed in terms of Fourier coefficients of the two-tangent form factors F£ b e , the coefficient reads 



*o (S) = 5(0) exp 



* l P k<A 



E E K^JeUe + 1) " W6»% fi 



v e of v 



For to be a well-defined state in the gauge-invariant Hilbert space Hq, the norm 



l*o|| 2 =AA 2 E I5(0)| 2 exp 



scr; 



E 

Zt P fc<A 



E E F v b e(k)je(je + 1) " V^%,0 



(107) 



(108) 



has to be finite. 

The number of possible spin network graphs on T£ is finite. Therefore, in order to show that (108) is 
finite, it suffices to prove that for every graph 7 on 7^, 

2" 

E E K b e(k)jeUe + 1) " WS a %, 



E IWPexp 



S with 
graph 7 



W2 E 



-u e of -u 



< OO . 



(109) 



Let us abbreviate j e (j e + 1) by c e . The factor 5(0) gives a polynomial in the j e 's that depends on the 
connectivity of the graph 7. For a suitable polynomial P[c e ] in c e , we have 

|5(0)| 2 <P[c e ]. (110) 



Thus, an upper bound on (109) is given by 

II E ) P \-°e] eX P ( - E C eCe' + E N e C e + Kj , 

e j e =l/2/ \ ee' e ) 

where M ee i is a symmetric positive matrix. By a linear transformation 
we render M ee i diagonal and arrive at the upper estimate 

IJ E I exp ( - E A eC 2 + E NeCe ) , 

e c e =0/ V e e / 



(111) 



(112) 



(113) 



where the eigenvalues A e are positive and P is a polynomial in c e . By using further estimates with integrals, 
(113) can be shown to be convergent. 

We conclude that the state (104) is an element in the gauge-invariant Hilbert space TLq. 
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3.5 Inclusion of phase factor 



Equation (104) is not yet the final form of the LQG state: what is missing is the contribution from the 
phase factor 

(114) 



exp 



^/d 3 x e^{x)e acd d c et{x) 



that we have ignored so far in our transition from Fock to loop state. The treatment of the additional 
integral is quite analogous to that of the Gaussian term. There is only one difference — a difficulty that 
arises when implementing the differential operator d c on the triangulation T£. If we were using a hypercubic 
lattice, we would simply get 



* {S) = 5(0) exp 



1 



d :i x / d*y W A (x - y) hf(x) hf(y) 



'ab/ 



(115) 



where V c stands for the lattice derivative in the c-direction. On a triangulation an implementation of 
d c is more complicated and requires additional weighting factors to take account of the geometry of the 
triangulation. We do not determine this in detail and content ourselves with saying that the final state has 
the wavefunction 



(5) = 5(0) exp 



1 



d 3 x / d 3 y W A (x - y) hf(x) hf(y) + phase term 



where the phase term is an analogue of 



— l — Y hf(v) e acd V c hf(i 
4HkP V 



(116) 



(117) 



on a triangular lattice. 



4 Graviton states 



The steps that led us from the vacuum of linearized gravity to the state can be repeated in complete 
analogy for gravitons. Let us go back to the Schrodinger representation of linearized extended ADM gravity 
(see 2.2). For k / 0, we define creation and annihilation operators 



ai(k) :-- 



^ etik) + 7^k 



k 



such that 



ai(k),a](k') 



Vhnk 



Ki{k) , 



The one-graviton state with polarization i and momentum k reads 



As usual, we can write this also with tensors: 

Vi,k[e2 d (k)]=2 




e?(k)eT(k)ne~ b d (k)} 



(118) 



(119) 



(120) 
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The canonical transformation to Ashtekar-Barbero variables adds the phase factor 



2i 



exp(^X>ei(£0e 2 (A0 



(121) 



Next we extend the functional from the reduced to the full configuration space, as we did in sec. 3.1. This 
gives us 



k 



*i,ki E m} = 2\lj-e?(k)eUk)*[E*(k)] 



2^e<?(k)-L J d 3 x e**e fa (s) ««(*)] 



(122) 



We replace ei a (x) by \ (x) — S la ^j (using again the argument that higher orders in e can be neglected), 
and arrive at 



(123) 



We bring this into the complexifier form, make the transition to Ti and finally apply the gauge projector. 
The result is the state 



x exp 



ATik 



d 6 x I d A y W A (x - y) hf(x) hf(y) + phase term 



(124) 



in the gauge-invariant Hilbert space Hq. We define an associated wavefunction 



% A (S) := S(0)\l-e?(k)hr(h) 



x exp 



AflK 



d 6 x I d 6 y W A (x - y) hf(x) hf(y) + phase term 



and write the state as 



(125) 



(126) 



SCT* 



In the same way, we construct multiply excited states. Denote the polarizations and momenta of the 
gravitons by i±, k_i, . . . ; iN,k.N- Then, the iV-graviton state in TLq becomes 



^'ii,fc 1 ;...;ijv,fc ]v — ■^'ii,k l ;...;i N ,k N E ^ii,k l ;...;i N ,k N (S)S ' 



(127) 



where 



A' 



k u 



,n=l 



x exp 



4hn 



d 6 x / d 6 y W\(x - y) hf(x) hf(y) + phase term 



(128) 



The normalization factor ftfi 1> k 1 -,...-,i N ,k N depends on the excitation number of each mode. 
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5 Semiclassical properties of the vacuum state 



By construction, the form of the vacuum is similar to that of the Fock space functional we started 
from. The exponential is still of a Gaussian type, where now the role of the fluctuation variable is played 
by the spin networks S and their associated (inverse densitized) metric gg b . As before, these fluctuations 
are non- locally correlated by the kernel W\. It is immediate from (116) that for most spin networks the 
coefficient ^o(S) is exponentially damped. An absolute value of the order 0.1 to 1 is only attained for a 
relatively small class of spin networks: by analogy with quantum mechanics, we say that these spin networks 
constitute the "peak region" of \I/o- The "position" of the peak itself is given by those spin networks for 
which ^>o(S) is exactly 1. 

The state inherits its Gaussian property from the semiclassical peakedness of the original Fock state. 
This suggests that we interpret the peakedness of as the way in which semiclassicality manifests itself on 
the level of spin networks: i.e. we interpret spin networks in the peak region as semiclassical fluctuations, 
and the spin networks at the peak position as the classical configuration. 

From that point of view, it would be interesting to know where exactly the peak is located; that is, for 
which spin networks S the coefficient \E f o(5') reaches its maximum. Below we analyze this question and try 
to estimate the peak position: at first for general values of the cutoff length, and then, in section 5.2, for 
the limit where l\ is much smaller than the Planck length. 



5.1 Peak position 

We have three length scales: the Planck length l p , the length cutoff l\ = ir/A corresponding to the momen- 
tum cutoff A, and the size of the 3-torus L = V 1 ^ 3 . The Fourier coefficient of the two-tangent form factor 
is 



*&Qs) = ^= (^j £ ds £ ds> e a (s)e b (s')S A (0) e"^ 



(129) 



The factor 5\ (0) gives 



<5a(0) 



V 



k<A 
1 4 



V (2n/L) 3 



d 3 k 



k<A 



-vrA 3 = -- 
(2vr)3 3 6 ll 



(130) 



In section 3.3, we assumed that the edges of the triangulation are straight and that their lengths are more 
or less equal to I a- Hence we can approximate (129) by 



F v ab e (k) 



Here, rj stands for 



v ■= —p 2 L~ 3 ll I7 l 
i 96 y v A 

and n e denotes the normalized direction vector of the edge e: 

e 



n,, 



(131) 
(132) 

(133) 



The j e -dependence of the S(0)-factor in (107) is polynomial. When determining the peak of the state, we 
can neglect it relative to the exponential dependence. Therefore, the peak condition becomes 

2 



fc<A 



Ce n a e n b e e-^ - 5 ab 5 kfl 



v e of v 
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V v e of v ) 0<k<A 

It is convenient to introduce the matrix 



2 

= minimal. (134) 



M ab :=Y, E VCeKnl (135) 

v e of v 

its traceless part 

T ab := M ab - - M 5 ab , (136) 
3 

and the function 

f ab (x):=WJ2 E ^ e <n»%-^). (137) 

v e of v 

With this notation, condition (134) takes the form 

, 2 



u (M ab -5 ab ) + ]T "(k)\f ab (k)\ 2 

0<k<A 

= 3lo (\m - l) 2 + u T ab T ab + J2 l/ a6 (^)! 2 = minimal. (138) 

^ 3 ' o<fc<A 

If we minimize each term separately, we obtain the three conditions 

(a) T ab = 0, i.e. M ab is isotropic. 

(b) M = ]T 2 Vc e = 3. 

v e of v 

(c) / a6 (A;) = for < k < A. 

Condition (a) requires that on average the spin network edges are isotropic; (b) requires a certain mean 
value for c e = j e (j e + 1)> an d (c) demands (weighted) homogeneity of spin network edges up to the scale 

Let us first consider (b): we define the mean value of c e on 7^ by 

s A = ^E c e > ( 139 ) 



e 



where iV e is the total number of edges of the dual complex Noting that N e equals twice the total 
number N v of vertices of 7^, we can write condition (b) as 

ca = • (140) 

When a triangulation consists of regular tetrahedrons of side length a, the dual edge length a* and the 
tetrahedron volume Vr are given by 

a* = ^ga, V T = y^a 3 = V3al. (141) 
The tetrahedrons of our triangulation are very close to being regular, therefore 



21 



25 
20 
15 

< 

10 
5 

°0 0.5 1 1.5 2 

/A 

Figure 1: Dependence of j\ on cutoff length l\. 

and 

3 96 2 LPl\ y/3Z^ 
« . (143) 

If we define the number j\ by 

jA(jA + l):=40/r 2 (^ , (144) 

condition (b) becomes 

CA=jA(jA + l). (145) 

For certain values of I a, the "spin" j\ is a half- integer (see Fig. 1). In these cases, it is easy to satisfy (a), 

(b) and (c): take the entire dual complex T£ as the spin network graph and label all edges with the spin j\. 
Clearly, the mean value of c e is ja(j'a + 1), so (b) is fulfilled. Since we have chosen a very homogenuous and 
isotropic triangulation, the uniform spin distribution over all edges of 7^ automatically implies isotropy (a) 
and homogeneity (c). 

Thus, for values of Ia which yield a half-integer j\, the peak of the state (116) lies at spin networks 
which have the entire dual complex as graph and all spin labels equal to j\. This peak exhibits a degeneracy 
due to the remaining freedom in choosing intertwiners at vertices. We see from (144) that the required spin 
depends strongly on the length cutoff I a- beyond I a > Q.7\f]3lp it grows quadratically in I a- In the range 
QAyJJilp < Ia < 0.7\fj3l p , it takes values between 1/2 and 10. 

What are the peaks of the state when j'a lies between two half- integers j\ < j-p. In that case, we 
have two opposing tendencies: condition (b) requires a mean value ca = jk{jh + 1) between + 1) and 
J2O2 + 1)) and therefore at least an alternation of labels between j± and ji- This is in conflict with condition 

(c) , which demands spatial homogeneity on all scales larger or equal to I a- That is, we cannot minimize all 
terms in (134) simultaneously. 

There are essentially two possibilities how the peak state could behave: one is that it becomes some 
blend-over of the uniform spin distributions for ]a = ji and ja = J2- The other possibility: the configuration 
keeps the uniform spin j\ while j\ moves away from j\ towards j2, and at some point it makes an abrupt 
transition to the uniform j'2-configuration. 
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5.2 Limit Z A -C l p 

Instead of discussing this in detail for general I a, we will concentrate on the most interesting case, namely, 
when I a is much smaller than the Planck length. Then, j\ lies between and 1/2. To determine the peak 
states, we introduce three classes of spin networks that we call Sb, Sp and S a b c - 

Sb denotes the class of all spin networks on that meet condition (b), whereas Sp stands for those 
which violate (b). S a bc consists of the spin networks which satisfy (a) and (b), and minimize the violation 
of condition (c), i.e. they minimize the third term in (138), while the first two terms are zero. 

We will first determine the properties of spin networks in S a b c , and then argue that this class of states 
is the minimizing solution for condition (138). 

We begin by analyzing condition (b). Suppose that for some A' < A, there is a triangulation of the type 
T\> whose dual complex T£, is coarser than T£ (for the definitons, see p. 14). Consider the restricted mean 
value 

(146) 



CA' 



where the sum runs only over edges of the coarser complex T£,, and N e > is the total number of its edges. By 
repeating the steps that led up to (145), we find that for spin networks on T£,, condition (b) is equivalent 
to 

c A ' =JA'Ua' + 1), (147) 

This equation generalizes equation (145) in the sense that it expresses (b) as a constraint on various scales of 
T^. In particular, for uniform spin distributions, it gives a relation between spin and scale of the distribution. 

Using this reformulation of condition (b), we immediately find a large class of spin networks that solve 
it: according to (147), a spin network meets condition (b), if it has a dual complex C T£ as its graph and 
a uniform spin labelling with ja>- For given ja>, we call the set of such spin networks Sj A ,. The admissible 
scales A' are those where j\i takes a half-integer value. For example, j\i = 1/2 corresponds to the scale 
0.4^/5^. By definition, a complex T£, is nearly isotropic, so spin networks in Sj A , meet also condition (a). 

The third condition on S a b c states that the spin networks should minimize the term 



Z = E ^)\f(k)\\ 

0<k<A 



(148) 



which measures inhomogeneity. A state in the class Sj A , breaks homogeneity (condition (c)) at the scale 
A'. Using (137) and (144), we see that it yields 



i = E <*>(*) 



0<fc<A 

A' 
A' 



E E iic e n a e n b e e *z 

v e of v 

jiWir^A'CM' + i) 




_A. 

Ia 1 



A' 



(149) 



In going from the third to the fourth line, we used that rjN v ~ j\ 1 (jA + l) -1 - Loosely speaking, high 
momentum wins over small momentum scale, and of all the sets Sj A , it is the one with the highest momentum 
scale A' that minimizes (148). Clearly, that set is S 1/2 - 
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Given the domination of low spin in (149), it appears unlikely that configurations with mixed spin labels 
lead to the same or a lower value. If that expectation is correct, S 1/2 does not only minimize the sequence 
S 1/2 , S 1 , 5*3/2, • • ., but also all other states that satisfy (a) and (b). That is, S 1/2 equals S a b c - (When saying 
so, we are slightly imprecise: the minimizing conditions do not require that the graph takes exactly the 
form of a dual complex of a triangulation. The set S a b c can also contain spin networks with other graphs 
Gh< as long as they are isotropic, homogeneous and have a length scale corresponding to j\i = 1/2. When 
speaking of S 1/2 in the following, we mean to include these additional spin networks.) 

Let us now consider the relation of S a bc to the entirety of configurations on T A *: we distinguish between 
configurations that satisfy (b) (the class Sb) and those which do not (the class Sp). By construction, S a bc 
is the subset of Sb that comes closest to the minimum. Regarding Sp, there is, in principle, the possibility 
that the violation of (b) is compensated by greater homogeneity of the spin distribution. To check this, we 
compare expression (138) for S a b c and Sp: for S a b c = S 1/2 we just get the inhomogeneity, i.e. 

Z-aJW 2 -^. (150) 



IaJ I 



On the side of Sp, we obtain 



^J l -M-l) +u>»T ab T ah + £ u(k)\f{kf 

Vd J 0<k<A 

> 3w Yl ^Ce ~ 1 ) 

V v e of w / 



/4 ^ 2 
~ LO I g V N v c A - 1 

~ ^° 1)2 ^ A ~ jA ^ jA + ^ ( 151 ^ 

~ ^o(^J [ca-JaUa + I)} 2 ■ (152) 

If we think about states in Sp whose mean value ca has a finite difference to ja(ja + 1) 5 the term (152) 
blows up much faster than (150) as l\ — > 0, so the S 1/2 states are preferred. Alternatively, we could consider 
states Sp that have a fixed relative difference between c\ and j'aC/a + 1), so that (151) remains constant. In 
that case, the inhomogeneity becomes the dominating criterion. Since j'a(ja + 1) gets close to zero, ca does 
so too, and the states Sp become similar to Sb states. In fact, there is no obvious reason why such states 
should have advantage over the S^-states. Hence we expect again that the S 1/2 class provides the lowest 
value. 

This would mean that for /a *C l p , the peak of the state consists of spin networks with the following 
properties: they are uniformly labelled by 1/2, their graphs lie on T£, are homogeneous and isotropic, and 
the effective length scale of the graphs is close to \fj3l p , regardless of how fine T£ is. 



6 Summary and discussion 
6.1 Summary of results 

The basic idea of our approach is to start from the free Fock vacuum of linearized gravity and construct 
from it a state that could play the role of the "free" vacuum in loop quantum gravity. In making the 
transition from Fock to loop state, we have to take various choices that relate the field variables of the 
former to the polymer-like degrees of freedom of the latter. We have done so with the aim of making the 
state reasonably simple while preserving, as far as possible, the physical properties of the original state. 
Let us repeat the logic of our construction: 
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1. We linearize extended ADM gravity around a flat background on T 3 x R. 

2. We apply a reduced phase space quantization and specify the free vacuum. 

3. We perform the canonical transformation to linearized Ashtekar-Barbero variables and implement it 
as a unitary transformation in the quantum theory. 

4. We regularize the state by a momentum cutoff A. 

5. The transition to the LQG Hilbert space is achieved in five steps: 

(a) We extend the state from a functional of reduced to a functional of full triads. 

(b) We replace the fluctuation ef = Ef - Sf by 

\ (g a \x) - 5 ab ) = \ (Ef(x)E!(x) - S ab ) , 

the fluctuation in the densitized inverse metric associated to Ef. 

(c) We bring the state into the complexifier form. 

(d) We replace the momentum cutoff A by a regular triangulation 7a of length scale n/A, and thus 
obtain a functional of generalized connections. 

(e) Gauge projection yields a state in the gauge-invariant Hilbert space Ho- 

Ideally, this procedure should be completed by an averaging over the 3d-diff group, so that one receives a 
state in TCdis- We have not carried out this step. 
The state we get is a superposition 

y =AT ]T Vo(S)S*. (153) 

SCT A * 

The sum ranges over all spin networks S whose graph lies on the dual complex T£ of the triangulation. The 
coefficients ^o(S) are given by 



*o(5) =5(0)exp 



— - / d 3 x / d 3 y W A (x - y) hf{x) hf(y) + phase term 
Ann J J — — 



(154) 



To each spin network S we associate an inverse metric, and h'g stands for the difference between g^ b and 
the flat background metric: 

hf = gf-S ab . (155) 
gg b is distributional and has only support on vertices of the spin network graph: 

~gf(x) = J2~gf(v)6(x-x v ) (156) 

V 

and 

~9s(v) = ^ lpA(0) £ j e (j e + 1) f 1 ds f 1 ds> e a (s)e b (s') . (157) 
lb eolv Jo Jo 

The functions W\{x — y) and 5a (x) are regularized forms of the kernel l/|x — y| 4 and the delta distribution 
respectively. On a hypercubic lattice, the phase term takes the form 

^Jh^)^^) (158) 
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where V c denotes the lattice derivative. Here, we are working with triangulations, so this formula has to be 
adapted to the geometry of the triangular lattice. In complete analogy, we have also constructed iV-graviton 
states. 

In sec. 5, we investigated the ^-dependence of the coefficient ^>o(S): it is of Gaussian type and expo- 
nentially damped for most spin networks. At the peak of the Gaussian, the labelling is determined by a 
characteristic "spin" j\ which depends strongly on the cutoff scale (see Fig. 1): 



We have shown that for cutoff lengths l\ where j\ is a half-integer, the peak spin networks have the entire 
dual complex 7^ as their graph and the edges are uniformly labelled by j\. For values of l\ where j\ lies 
between two nonzero half- integers ji < j2, the graph is again the entire complex and we expect that the 
spin labelling is either uniformly j±, j2, or a mixture of both. When l\ is much smaller than the Planck 
length, our analysis supports the idea that the peak spin networks become independent of the cutoff, have 
spin label 1/2 and homogenous isotropic graphs at a length scale close to \fj3l p . 

6.2 Relation to other approaches 

Let us point out some similarities and differences to what has already appeared in the literature: 

The idea of using linearized states is not new and has been investigated before by Ashtekar, Rovelli 
k, Smolin [13] and Varadarajan [14]. In both cases, however, the Fock states were transferred to a loop 
representation of linearized gravity, and not to the full representation, as we do here. Our state is also 
related to the class of complexifier coherent states which are used by Thiemann and collaborators [15]-[22]: 
what distinguishes our proposal is the particular choice of "complexifier" and the fact that we leave the 
background field outside the delta functional — a circumstance which facilitates the investigation of peak 
properties considerably. One could say that we keep the background in the "complexifier" , though in the 
strict sense of the word this object is no complexifier anymore. Thiemann also introduces cutoff states that 
have support on a fixed graph, but it is not clear which graph should be chosen or how such a choice should 
be justified. In our case, the cutoff graph is fixed on physical grounds, namely, by the requirement that it 
emulates the momentum cutoff of the original Fock state. 

Ashtekar and Lewandowski [23] have proposed a coherent state that is similar to ours in that it is 
based on J-insertions at vertices (see eqn. (95)). The way in which the J-operators (and form factors) are 
contracted is different, however. Moreover, their convolution kernel is l/k (up to regulators), while ours 
corresponds to k. This difference can be explained as follows: Ashtekar and Lewandowski start out from 
£7(1) Maxwell theory, where the Hamiltonian is H ~ E 2 + k 2 B 2 , and then generalize from U(l) to SU (2). In 
the E'-representation, the vacuum wavefunctional contains the kernel l/k and that is also the kernel which 
appears in the complexifier. Our state, on the other hand, comes from the linearization of gravity where 
the linearized Hamiltonian reads H ~ K 2 + k 2 e 2 . Thus, the vacuum functional in e has k as its kernel and 
inherits it to the loop state. 

Our results on the peak of the state clearly show a connection to the early weave approach [24]-[27], 
where a smeared metric operator was used to define eigenstates of the 3-metric. The construction of \I>o 
involves a smeared inverse metric operator 



and one of the conditions for maximizing ^o^) consists in the requirement that S is an eigenstate of (160) 
with eigenvalue equal to 5 ab . Thus, the peak spin networks of can be considered as a certain form of 
weave. For l\ > \fj3l p , these weave states have the entire dual complex as graph, and the cutoff-dependence 
of the spin labelling follows the same logic as in the weave approach: for a large cutoff length l\, the edges 
of the dual complex are dispersed far apart and they have to be labelled with large spins to generate a field 




(159) 




(160) 
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strength equivalent to the background. For smaller I a, the available edges are more densely distributed and 
a smaller mean spin suffices to attain the same field strength. It is somewhere near l\ = \f]il p where the 
mean spin reaches values of the order 1. 

If I a is considerably smaller than the Planck length, one would expect that the average j e has to be much 
smaller than the minimum nonzero spin 1/2, so that only now and then an edge carries a spin 1/2, while 
most other edges are labelled trivially. In other words, one could choose the triangulation increasingly finer 
than the Planck length, but the spin networks would stay apart as if they were trying to keep a triangulation 
fineness at a larger length scale. This is, in fact, what our analysis indicates: namely, that for l\ <C l p , the 
peak spin networks are in the fundamental representation and have isotropic homogenuous graphs whose 
edges maintain a distance scale close to the Planck scale. This stands in correspondence to the key result 
on weaves [24], which states that the approximation of 3-metrics cannot be improved by making the lattice 
finer than l p . Let us mention that in [28], Corichi and Reyes have also constructed a state that is peaked 
around weave-like spin networks in the fundamental representation. 

6.3 Continuum limit of free vacuum? 

The cutoff-independence of the peak suggests that the entire state could become independent of ^a when 
l\ goes to zero. Note that the limit should be taken as 

lim A/a Y ^ {S)S* (161) 

ScT* 

where A/a normalizes the state for each value ^a in the sequence. Since gg b diverges for Ia — ► 0, any deviation 
from the peak will be infinitely suppressed: the limiting state loses its spreading and becomes effectively a 
weave state. 

The reason for this behaviour lies in the fact the we have not replaced the momentum cutoff completely 
when going from the Fock to the loop state. Instead we ended up with a mixed scheme where part of the 
regularization is provided by the cutoff graph T£ and part of it by the momentum cutoff in the kernel Wa 
and the factor (5a (0) in (157). Such a mixed regularization may be ok for I a > l p , but it becomes physically 
questionable when the regulator is removed: then, the triangulation part leads to an effective cutoff at the 
Planck scale ("space is discrete"), while the momentum part lets 5a(0) diverge ("space is continuous"). 

To arrive at a more meaningful continuum limit, it would be desirable to express the regularization 
entirely in terms of the triangulation. Carlo Rovelli has suggested that one might achieve this by starting 
from a modification of the Fock state (83): observe that the kernel Wa and the delta function 5a require 
the background triad E c i as an input, so we write them as Wa,e c1 and 8a,e c1 - In the state functional, we 
replace these quantities by Wa,e and 5a,e, since it does not change the part of the exponential which is 
quadratic in the fluctuation E — E c \. The kernel Wa,e is the same as a/— Aa,e where Aa,e denotes the 
A-regularized Laplace operator associated to the triad E. 

The next step is the transition to the LQG state, that is, one has to decide on how to implement 
the propagator Wa,e = \/—Aa,e in the LQG complexifier. The first thing to note is that when making 
this transition, the A in Aa,_b can be dropped because the triangulation takes care of the momentum 
regularization. Therefore, what remains to be done is to find a loop quantization A^ of A^ on this 
triangulation. If things are simple, it will act diagonally on spin networks and contributes a factor W e e >(S) 
for every pair e, e' of edges of the spin network S on which the complexifier acts. In a similar way, one 
might be able to define a loop analogue of 5a (0). 

The conjecture is that by this procedure one arrives at a "free" vacum that becomes independent of the 
cutoff when I a goes to zero and keeps the form that it had for I a ~ l p . 

6.4 Free theory, perturbation theory and renormalization? 

We finally come to the big question that motivates this work: can one obtain a semiclassical perturbation 
series for loop quantum gravity, and if yes, how? 
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To start with, let us recall how perturbation theory works in ordinary QFT, or even simpler, in quantum 
mechanics. Consider the anharmonic oscillator with the Hamiltonian 



H = ^P 2 + ^ 2 (Q-Q ) 2 + X(Q-Qo) 4 - (162) 

Assume that we want to do perturbation theory around the static solution Q(t) = Qo, P(t) = 0. For that 
purpose, we introduce the relative variables 

q-=Q-Qo, p:=P, (163) 
and "expand" the Hamiltonian in q and p: 

H = ^p 2 + l -u J 2 q 2 + \ q i . (164) 

Perturbation theory rests on the idea that we consider only processes that involve a certain subset of states; 
namely those whose wavefunction is peaked near the classical position and momentum, and strongly damped 
farther away from it. If the width of these states is sufficiently small, one may ignore higher order terms to 
first approximation, and use instead the free Hamiltonian 

H = \v 2 + \^ 2 q 2 . (165) 

The associated free vacuum is 

MQ) =AAex p("|^ 2 ) • (166) 

For this approximation to be consistent, the width of should be small enough that the neglection of 
higher orders is admissible. We check this by applying the full potential V to ipQ-. 

V^( q )=M(^L0 2 q 2 + X q Ae W (-^qA . (167) 



The width of ipo is y/h/u, so the consistency condition becomes 



or 

u 3 



-J 2 - > A - , (168) 



*«-. (169) 

If it is satisfied, there is hope that a perturbative treatment leads to meaningful results. Otherwise we are 
in the non-perturbative regime. 

Let us transfer this logic to the linearization of ADM gravity on T 3 x R. When we keep all orders in 
the reduced fluctuation variables, the Hamiltonian reads 

H = -J d 3 x [(K l a K> - KiK l b ) (5 + e<j) («5f + ef) - |det (<5« + e<j) \ R + efj\ . (170) 

Suppose we quantize this system as we did in the free case, i.e. by a standard canonical quantization and 
using a momentum cutoff A as a regulator. Then, a calculation analogous to (167) yields the expansion 
parameter 

a: = Jf2 • ( 171 ) 
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Here, we have taken k k to be the gravitational coupling at the cutoff scale. L is the size of the torus and 
the time scale of the process we consider. The simplest way to relate k a to the classical value is dimensional 
scaling, i.e. 

I 2 

k a = ^k. (172) 



In that case, the parameter (171) becomes 



1 T 1 / 2 

a = ^ . (173) 

'A 



This estimate indicates that for L of the order 1cm, perturbation theory requires a cutoff length larger than 
10 -22 cm. For l\ = lp, the theory would be highly non-perturbative. 

Of course, what we actually want to analyze is loop quantum gravity, and not a conventional quantization 
of gravity with an arbitrary cutoff. With that aim in mind, we transferred the free vacuum of (170) to the 
LQG Hilbert space. The semiclassical spreading around e = turned into a Gaussian spreading around 
weave-like spin networks. The more the spin labelling of a spin network differs from that of the weave, the 
more it is damped in the sum (153). The quantity that measures the deviation from the peak is 

hg = gf-5 a \ (174) 

and stands in obvious correspondence to the fluctuation e of the traditional approach. Therefore, it appears 
natural to ask the following questions: could this correspondence give us a hint on how semiclassicality 
manifests itself in full loop quantum gravity? Could we take the size of /Y| b , or a similar quantity, as the 
loop analogue of the fluctuation, and use it as a measure for telling whether a spin network is a "small" 
or a "large" fluctuation w.r.t. a semiclassical state? Suppose that we let the full Hamiltonian constraint 
C act on a state like that is peaked around weaves. Would there be parts of C that are more relevant 
than others, in the same way that uj 2 q 2 /2 is more relevant than Xq 4 when acting with the potential on 
exp(— ujq 2 /2h)? In other words, can we write C as a sum 

C = C + Ci (175) 

such that to first approximation Ci^o can be neglected relative to Co^o? To answer that question, one 
would have to analyze if the quantity 

^o(S)CS (176) 

can be approximated for spin networks close to the peak weaves. I.e. for those spin networks whose spin 
assignments differ from the weaves on only few edges and only by little spin. 

Thus, we are led to an ansatz that appears very similar to the one used by Smolin when he analyzes 
string perturbations of causal spin networks [36]: in that case, the "differing" edges are chosen to form 
loops and are specified on an entire history of causal spin networks. The sequence of loop-like deviations is 
viewed as a string worldsheet, and the resulting variation in the amplitude is to some extent evaluated. 

The inherent cutoff-property of LQG suggests that loop quantization is related to conventional quanti- 
zation schemes where a cutoff l\ = lp is put in by hand. Similarly as there, the expansion parameter could 
be large and perturbation theory impossible. In that case, a coarse-graining procedure [34, 35] may be nec- 
essary to compute an effective Hamiltonian constraint for a lower scale l\, where the expansion parameter 
is small. Our state *o might be useful in such attempts because it is an element in the full Hilbert space 
and represents at the same time a semiclassical state at the cutoff scale I a- One could try to extract an 
effective matrix element 

(^ol^Al^o) from (* |P|*o), (177) 

where P is the "bare" projector onto the physical Hilbert space, and Pa is based on a renormalized Hamil- 
tonian constraint at the cutoff scale l\ < l p . 

Admittedly, the considerations of this last section are vague and need to be concretized in many ways. 
We hope to do so in future work. 
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